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Abstract 

State-selected Coulomb-crystallized molecular ions were employed for the first time in ion-molecule reaction studies 
using the prototypical charge-transfer process + N2 —> N2 + as an example. By preparing the react ant ions in a 
well-defined rovibrational state and localizing them in space by sympathetic cooling to millikelvin temperatures in an ion 
trap, state- and energy-controlled reaction experiments with sensitivities on the level of single ions were performed. The 
experimental results were interpreted with quasi-classical trajectory simulations on a six-dimensional potential-energy 
surface which provided detailed insight into translation-to-rotation energy transfer occurring during charge transfer 
between N2 and Nj. 
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1. Introduction 

Achieving full control over the internal quantum states 
as well as the collision energy of the reaction partners to 
probe the detailed state-to-state dynamics and energy de- 
pendence of chemical processes has been one of the long- 
standing aspirations in gas-phase chemical studies. In the 
past, a range of experimental methods has been developed 
to address these objectives, e.g., internal cooling of the re- 
act ants in supersonic expansions, experiments with crossed 
molecular beams at variable relative velocities and laser 
preparation of the reactants in specific quantum states 
(see, e.g., Refs. [1, 2, 3 and literature cited therein). Very 
recently, new techniques have been established which al- 
low the preparation of molecules at extremely low transla- 
tional temperatures < 1 K and at the same time enable an 
unprecedented degree of control over their kinetic energy 
[H|5j[6]. In conjunction with the simultaneous preparation 
of the molecules in well-defined rotational- vibrational and 
even hyperfine states, "cold molecules" methods are now 
starting to pave the way for studies of collisional processes 
and chemical reactions in new physical regimes and at lev- 
els of detail and sensitivity which have not been possible 
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Whereas most efforts have so far concentrated on neutral 

molecules, techniques for cooling and controlling molecu- 
lar ions have recently made impressive progress as well. 
From an experimental perspective, the sensitivity of the 
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ion motion to weak stray electric fields in the apparatus 
renders the precise control of their kinetic energy difficult. 
This problem can be overcome by trapping the ions and 
cooling them sympathetically to millikelvin temperatures 
by the interaction with co-trapped laser-cooled atomic ions 
[T2] . Under these conditions, the ions localize in space to 
form ordered structures usually referred to as "Coulomb 
crystals" in which it is possible to observe, manipulate and 
address single ions [6j [13] . 

Moreover, the preparation of molecular ions in well- 
defined quantum states is experimentally challenging and 
indeed only a handful of reaction studies with rotationally 
state-selected ions have been reported thus far (see, e.g., 
[HI HSJ HH [TTj HE]). Only very recently it has become 
possible to prepare Coulomb-crystallized molecular ions 
in well-defined internal quantum states [I9j [20j [2lJ [22] , so 
that a simultaneous control over the kinetic energies, posi- 
tions and internal states of the ions can now be achieved. 

In the present article, we report the first study of 
a chemical reaction with rovibrationally state-selected 
Coulomb-crystallized molecular ions. We investi- 
gate the prototypical symmetric charge-transfer (CT) 
reaction + N2 —> N2 + between state-prepared 
sympathetically-cooled Nj! - ions and internally cold N2 
molecules from a supersonic expansion at a well-defined 
collision energy. 

The mechanism and kinetics of this reaction have been 
studied previously, see, e.g., Refs. [23l f24| f25 | f26| f27 j [28] 
and references therein. At low energies, the reaction 
was found to proceed via a long-lived N|~ reaction com- 
plex which forms at the collision (Langevin) rate [24j 
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Figure 1: Schematic of the experimental setup. 

[25] [28] . Using isotopically labeled nitrogen ions, Frost 
et al. showed that near-thermal collisions between 
ions and N2 molecules in their vibrational ground states 
lead to a symmetric CT with a rate constant of k = 
4.2 x 10 -10 cm 3 s _1 which amounts to one half of the 
collision rate constant [25] . This observation was rational- 
ized in terms of a symmetric sharing of the charge between 
the two N2 moieties in the N|~ reaction complex so that 
the probability for CT amounts to 50% upon its breakup. 
Vibration-to- vibration (V-V) and vibration-to-translation 
(V-T) energy transfer occurring during CT have been in- 
vestigated in this system both experimentally and the- 
oretically, see Refs. 1 25] [26] EH 123 ED] and references 
therein. Kato et al. observed multi-quantum vibrational- 
energy transfer at near-thermal collision energies which 
was rationalized in terms of energy redistribution within 
the strongly bound N4 reaction complex [26]. Also, ev- 
idence for vibration-to-rotation and translation (V-R,T) 
energy transfer was found. 

To our knowledge, the present work is the first to achieve 
rotational state selection of the reactant ions and rota- 
tional resolution in the analysis of the product-ion state 
populations allowing to gain insight into the dynamics 
of translation-to-rotation (T-R) energy transfer occurring 
during CT. The Coulomb-crystal technology employed in 
the present study enabled us to monitor reactive collisions 
in small ensembles of 20 to 30 state-selected, spatially lo- 
calized ions with single-particle sensitivity. The experi- 
mental results were interpreted with quasi-classical trajec- 
tory calculations on a fully-dimensional potential-energy 
surface (PES) which provided insight into the charge- and 
energy-transfer dynamics underlying this fundamental re- 
action. 

2. Experimental methods 

The experimental setup for generating quantum-state 
selected Coulomb-crystallized ions (see Fig. [TJ has 
been described in detail in previous publications [19] 
[20]. Briefly, N+ ions in the rovibronic ground state 



X + 2 S+,v + = 0,7V + = 0, Fi were prepared by a 
[2 + V] resonance-enhanced threshold-photoionization se- 
quence via the a" 1 E+,v / = 0, J 7 = 2 intermediate level 
of neutral N2 as shown in Fig. [2] (a). Here, v + and v r 
(7V + and J') stand for the vibrational (rotational) quan- 
tum numbers of the cationic ground and neutral intermedi- 
ate states, respectively, and F\ denotes the spin-rotational 
component. The ions were generated inside a linear 
quadrupole ion trap [6[ [13] from a collimated, pulsed su- 
personic molecular beam of pure N2 gas. The molecu- 
lar beam was positioned at a distance of « 400 /im from 
the Coulomb crystal in order to prevent collisions between 
ions and neutrals during the loading phase. The rota- 
tional temperature of the N2 molecules in the beam was 
determined to be T rot « 10 K by resonance-enhanced 
multiphoton-ionization (REMPI) spectroscopy [20 corre- 
sponding to populations of 50, 25, 22, and 3 % in the 
J = 0,1,2,3 rotational states. The beam velocity was 
estimated to be ~ 787 ms _1 from flow-dynamics models 

ED- 

Immediately after their generation, typically 25 state- 
selected ions were sympathetically cooled to transla- 
tional temperatures of « 10 mK by the interaction with 
laser-cooled Ca + ions to form bi-component Coulomb crys- 
tals [6] [13] . Ca + ions were produced in the center of the 
trap by non-resonant photoionization of Ca atoms ema- 
nating from a Ca oven. The Ca + ions were laser cooled on 
the As 2 S\ 1 2 — >• 4p 2 Pi/2 transition using diode-laser radi- 
ation at 397 nm. Another diode-laser beam at 866 nm was 
used to repump population on the 3d 2 D%/2 4p 2 Pi/2 
transition to close the laser cooling cycle. The resulting 
Ca + /N2" bi-component crystals were imaged by collecting 
the spatially resolved laser-cooling fluorescence of the Ca + 
ions using a microscope coupled to camera. In the images, 
the positions of the non-fluorescing molecular ions are vis- 
ible as a dark region in the center of the crystals, see Fig. 
1(a)- 

Reactive collisions between state-selected ions and 
neutral N2 molecules were initiated by overlapping the 
molecular beam with the bi-component Coulomb crystal 
immediately after ion loading and sympathetic cooling. 
The collision energy E co \ ^0.045 eV in the present ex- 
periments was entirely dominated by the kinetic energy 
of the N2 molecules in the beam. Neglecting the small 
initial rotational excitation of neutral N2, the maximum 
kinetic energy of the product ions cannot exceed the to- 
tal kinetic energy of the reactants which is much smaller 
than the trap depth (>2 eV). Consequently, the product 
ions remained trapped and were sympathetically re-cooled 
into the Coulomb crystal. Because the product ions were 
chemically identical to the reactant ions, their presence 
manifested itself in a time-dependent increase of the pop- 
ulation in rotationally excited states in the ensemble of 
sympathetically-cooled ions. Vibrational excitation of 
the products was precluded on energetic grounds. Rate 
constants were determined by measuring the Nj! - spin- 
rotational state populations as a function of the effective 
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Figure 2: (a) Resonance-enhanced threshold-photoionization scheme 



used to generate Njj~ ions in the X+ 



0,AT+ = o,Fi 



spin-rovibrational ground state, (b) Laser-induced charge-transfer 
(LICT) scheme for probing the populations in the spin-rotational 
levels of the vibronic ground state of produced by CT 

reactions. See text for details. 



time of reaction with N2 molecules from the beam and 
fitting the results to a kinetic model as detailed in Sec. [4j 
The rotational-state populations of the product 
ions were probed by optical pumping to vibrationally 
excited levels to promote CT reactions with Ar atoms 
which are energetically forbidden in the vibrational ground 
state (laser-induced charge-transfer (LICT) spectroscopy 
pjJJ [20] ) . LICT was initiated by exciting transitions to 
selected rotational levels of the X + 2 S+ v + = — » 



A+ 2 IL, 



4 state with a pulsed dye laser operating 



at 613 nm, see Fig. [2] (b). Subsequent fluorescent de- 
cay lead to the population of vibrationally excited states 
v + > 1 in the electronic ground state with a probability 
of 93%. The removal of vibrationally excited Nj! - ions by 
CT with Ar gas leaked into the chamber was directly ob- 
served by the reduction of the core of non-fluorescing ions 
in the images, see Fig. H (b). The number of N+ ions lost 
by LICT represented a direct measure of the population 
in the initially excited spin-rotational state. This num- 
ber was determined by a comparison of the experimental 
fluorescence images with simulated images generated from 
molecular-dynamics simulations [SHU], see Fig. [3] (a). 

3. Theoretical methods 

3.1. Potential energy surfaces 

The ground state six-dimensional potential-energy sur- 
face (PES) of the N 2 -N+ system was computed at the 
UCCSD level [32j [33] with Dunning's correlation consis- 
tent polarized valence basis set (cc-pVTZ) [34 using Gaus- 
sian03 [35]. Single-point calculations were performed for 
5565 non-equivalent geometries with energies up to 1.9 eV 
relative to the minimum of the complex on a non- 
equidistant rectangular grid including the N-N separations 
7*1 and r2 in N2 and Nj, respectively, the center of mass 
distance R between the two diatomic molecules, and the 



angles #1, 62 and (f) (see the inset in Fig. |6|. The calcu- 
lated well-depth (D e ) of the complex is 1.254 eV, in good 
agreement with previous even higher-level theoretical es- 
timates (1.26 eV at the (RCCSD(T)/vqz(spdfg) level of 
theory [36]). For comparison, the experimentally deter- 
mined dissociation enerev of N+ is D = 1.06 eV [37]. 

The global PES was represented by three surfaces, re- 
ferred to as one "bound" (N4") and two "unbound" ones 
(N2-NJ and NJ-N2). They were distinguished by the dis- 
tance R and the localisation of the majority of the charge. 
For R < 7.09 ao (3.75 A) the system is considered to be 
bound whereas for R > 7.56 ao (4.0 A) it is unbound. The 
intermediate region (7.09 — 7.56 ao) represents the tran- 
sition region during the dynamics, where the surfaces are 
connected by a smooth switching- function [38] . 

To allow a) bond formation between any two atoms of 
N2 and and b) dissociation of the complex into either 
the charge-preserving or charge-transferring state at least 
8 force fields (FFs) are needed for the bound state and 4 
FFs are necessary for each of the unbound states. These 
FFs are related to each other through permutation of their 
parameters. Within each set of FFs always the lowest 
energy surface is followed, except for when they are close in 
energy (within ^0.02 eV) in which case they are smoothly 
joined by an energy difference-based switching function 
[38]. 

Morse and Lennard- Jones potentials were used for the 
bond stretches and the intermolecular interactions, respec- 
tively, for the unbound states. Morse parameters for the 
isolated molecules were determined from PES scans at the 
UCCSD level of theory. For the bound state, Morse poten- 
tials were used for all bonds of the complex and Lennard- 
Jones and electrostatic potentials between the edge atoms. 

The force field parameters were fitted to the 5565 ab 
initio points with a simplex algorithm [39 . The ab ini- 
tio energies for the unbound region were reproduced with 
a root-mean-square-deviation (RMSD) of « 0.013 eV. It 
was considerably more difficult to obtain a good fit for the 
bound state due to strong angular dependence of the po- 
tential. The final fit had an RMSD of w 0.06 eV by using 
two minimal symmetrized sets of FFs (twice 8, altogether 
16 FFs). Figure [6] shows an example of the quality of the 
fit. The well depth of the bound state was exactly repro- 
duced by the global surface. In the following, the bound 
state FFs are numbered 1-16, whereas those for the un- 
bound states are labeled 17-20 and 21-24, respectively. 

3.2. Molecular- dynamics (MD) simulations 

For the quasi-classical trajectory calculations a code 
with provisions for adiabatic reactive MD (ARMD) 
[40] was used. The Hamiltonian equations of motion 
were solved in Cartesian coordinates using the adaptive 
timestep Modified Extended Backward Differentiation 
Formulas method (MEBDFSO) [41]. 

The dynamics was initiated in the unbound state (FFs 21- 
24). Between centre-of-mass separations of R = 7.56 ao 



3 



and R = 7.09 ao the momentarily active unbound surface 
was smoothly switched to the bound FFs (numbers 1-16). 
The complex is considered to be formed for R < 7.09 ao 
and the dynamics is continued in the bound state. When 
a separation of R = 7.56 ao is reached again, it is de- 
termined which unbound FF is lowest in energy at the 
corresponding geometry and the dynamics is followed on 
this FF by smoothly switching from the bound state. Fol- 
lowing the dynamics further, the system either can return 
to the bound state (recrossing) or it can decay into ei- 
ther N 2 +Nj (charge preserving, FFs 21 to 24) or N^+N 2 
(charge transferring, FFs 17-20). The dynamics is followed 
until either the lifetime of the bound state exceeded 100 ps 
or the initial separation of the two fragments was reached 
again. Once the complex survives a few vibrational pe- 
riods, energy is expected to have randomized completely. 
Our results show that if the complex survives for at least 
0.4 ps the correlation between the reactant and product 
states is negligible. 

3.3. Analysis of final sates 

As a consequence of using classical dynamics, prod- 
ucts are formed with rovibrational energies and angular 
momenta which correspond to fractional quantum num- 
bers. Also, products having less than zero-point vibra- 
tional energy (ZPE) can be formed. This ZPE leakage is 
a shortcoming of quasi-classical simulations and various 
methods were proposed for either avoiding or correcting it 
[42j [43j 04] . In the present study, only trajectories which 
satisfy certain constraints in the product states, i.e., the 
total ZPE of the two molecules should be conserved within 
±10%, were further analyzed. 

4. Results and Discussion 

4.1. Product rotational- state distributions 

Fig. [3] (b) shows LICT measurements probing the pop- 
ulations in the 7V + = 0, Fi and = 1, F\ spin-rotational 
states as a function of the effective time of reaction with 
neutral N2 molecules from the molecular beam. For the 
N + = state, a marked decrease of the LICT efficiency 
can be observed with increasing reaction time, indicating 
the removal of the initially prepared N^, = ions by 
CT with N2. Conversely, an increase in the LICT efficiency 
out of the N+ = 1 state can be observed, indicating the 
generation of ions in rotationally excited states as a con- 
sequence of the reactive collisions. 

LICT measurements on the population of the = 
0,1,2 Fi 5 2 levels as a function of the reaction time are 
shown in Fig. [4] Two independent LICT measurements 
were performed for each spin-rotational level at four differ- 
ent effective reaction times. As a general trend, a decrease 
of the population in Af + = in favor of an increase of the 
population in A/" + > 1 was observed. This result indicates 
that product ions over the whole range of rotationally ex- 
cited states probed in the present study were generated 
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Figure 3: (a) Experimental false-color fluorescence image and 
molecular-dynamics simulation of a Ca + /N^~ bicomponent Coulomb 
crystal before reaction. The spatial distribution of the non- 
fluorescing molecular ions has been made visible in green in the 
simulated image, (b) LICT experiments probing the populations 
in the N+ 7V+ = 0,Fi and 7V+ = l,F± spin-rotational states. The 
panels show fluorescence images obtained after LICT as a function 
of the effective time of reaction t with neutral N2 molecules from the 
molecular beam. The LICT efficiency (corresponding to the number 
of non-fluorescing ions removed from the crystal in the boxed areas) 
decreases in = and increases in N~*~ = 1 with increasing reac- 
tion time, indicating the generation of rotationally excited Njj" ions 
from the CT reaction NJ+N2. See text for details. 

as a consequence of CT between NJ+N2. Because the 
neutral products leave the trap, no conclusions on their 
rotational-state distribution can be drawn. 

The observed redistribution of population in the ensem- 
ble of Coulomb-crystallized ions is the result of a se- 
quence of CT reactions with neutral molecules from the 
beam. The initially state-selected ions undergo CT 
collisions with rotationally cold N2 molecules. The result- 
ing product ions are sympathetically cooled into the 
Coulomb crystal whereupon they in turn can undergo colli- 
sions. In this way, the ensemble of originally state-selected 
N2" ions in the Coulomb crystal is lost and replaced with 
rotationally excited product ions. 

Close inspection of the data in Fig. [4] reveals two im- 
portant details. First, already in the beginning of the ex- 
periment when most of the reactive collisions occur with 
ions in iV + = 0, the generation of ions in all rotational 
levels probed in the present study can be observed. This 
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Figure 4: LICT efficiencies reflecting the populations in the lowest 
five spin-rotational levels of N^" as a function of the reaction time 
with neutral N2 molecules. Over the course of the reaction, the 
initial ensemble of state-selected ions is replaced by rotationally 
excited product ions. The lines are a fit of the data to the kinetic 
model shown in Fig. [5] (a). 

observation suggests that the production of product ions 
is feasible over a broad range of rotational states by reac- 
tions with in iV + = 0. The relevant rates appear to be 
of a similar magnitude, at least over the range of product 
states probed in the present study. Second, for a specific 
rotational level 7V + , there appears to be a slight preference 
for the generation of the F\ spin-rotational component in 
comparison to the F2 component. 

4.2. Kinetics 

The data were analyzed in terms of kinetic models tak- 
ing into account CT collisions of the initially prepared 
N2" (iV + = 0) reactant ions as well as consecutive reac- 
tions of the product ions, see Fig. [5] Because the 
number density of neutral N2 in the molecular beam is 
much larger than the number density of the ions in 
the trap, pseudo-first-order kinetics was assumed in both 
models. 

Based on the observation that all product rota- 
tional states probed in the present study (iV + < 2) 
seem to be produced with similar rates, it was as- 
sumed in a first model that all processes of the form 
N+(7V+)+N 2 (J) N 2 (J)+N+(7V+) (bars indicate the 
rotational quantum numbers after the decay of the reac- 
tion complex) occur with the same /c, see Fig. [5] (a). The 
complete loss of population from the manifold of states 
with 7V + < 2 was taken into account by an effective loss 
rate constant k' . In the absence of detailed information on 
the neutral product state distribution, k and k' represent 
effective rate constants averaged over all neutral product 
states. The spin-rotational level populations Ui were ob- 



Figure 5: Schematic of kinetic models used to fit the redistribution 
of rotational populations in the ensemble of sympathetically-cooled 
N^" ions as a consequence of CT collisions. Model assuming (a) an 
uniform rate constant k and (b) different rate constants k\ and k,2 
for processes connecting unlike and like spin-rotational levels F\p 
in the reactant and product ions, respectively. See text for details. 
gi denotes the degeneracies of the spin-rotational levels including 
nuclear-spin statistical weights. 

tained from the set of rate equations 

~^jT = k ^2(9i n j ~ 9jni) - k'm, (1) 

where gij represent degeneracy factors (see Fig. [5|. The 
initial state populations rti(t = 0) were also treated as a 
fit parameter, accounting for the uncertainty in determin- 
ing the starting time of the measurement following the re- 
alignment of the molecular beam after ion loading (see Sec. 
[2]). The fit of this model to the experimental data yielded 
the pseudo-first-order rate constants k = 0.0023(4) s _1 
and k' = 0.04(3) s _1 . Because the density of neutral 
molecules in the molecular beam is not precisely known, 
it was not possible to obtain second-order rate constants 
from these results. Previous rate measurements with ther- 
mal samples of N2 gas leaked into the chamber [20 were 
consistent with the result obtained by Frost et al. [25] that 
the total CT rate amounts to one half of the Langevin col- 
lision rate k L = 8.3 x 10 -10 cm 3 s _1 . 

The time-dependent populations of the iV + = 
0, l,2,Fi 5 2 states computed with the kinetic model using 
the fitted parameters are shown in Fig. [4] The agreement 
between calculated and experimental level populations is 
satisfactory given the spread of the measured level pop- 
ulations. The agreement vindicates the assumption that 
the state-specific rate constants do not (or only weakly) 
depend on the rotational state of the product ion within 
the uncertainty limits of the present measurement and the 
range of states studied. The kinetic model also reproduces 
the preferential generation of ions in the F\ compared to 
the F2 spin-rotational components of the same 7V + state 
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Figure 6: Comparison of the UCCSD/cc-pVTZ (black) and fitted 
energies (red) for linear structures with r\ = T2 = 2.08 ao as a 
function of R. The inset shows the relevant coordinates scanned in 
the ab initio calculations. 



suggesting that this effect is caused by the higher statisti- 
cal weight of the F\ components. 

To check the validity of this model, its main simplifica- 
tion, i.e., the restriction to a "universal" st at e-to- state rate 
constant /c, was relaxed. In a second kinetic model (see 
Fig. [5](b)), different rate constants k± and k^ for reactions 
connecting states with different and like spin-rotational la- 
bels Fi 5 2, respectively, were assumed to account for the 
preferential production of product ions in the F\ levels. 
The fit yielded rate coefficients k\ = 0.0023(9) s _1 ,/c2 = 
0.0022(7) s -1 and k' = 0.04(3) s _1 . The results obtained 
for ki and &2 agree with each other and with the value of 
k within the uncertainty limits, supporting the conclusion 
that the increased production rates of ions in the F\ lev- 
els is indeed solely caused by the higher statistical weights 
associated with these channels. 

4-3. Trajectory calculations 

Quasi-classical trajectory calculations were performed 
to interpret the experimental findings and gain further 
insight into the energy-transfer mechanisms involved in 
the reaction. Initial states for the trajectory calculations 
correspond to semi-classically quantized states of rotating 
Morse oscillators [45] . In the experiment the molecules are 
exclusively in their ground vibrational state (v = 0). The 
initial rotational angular momentum 7V + of ions was 
set to zero as prepared in the experiment whereas that of 
N2 was randomly sampled according to the experimentally 
measured distribution (see Sect. Initial coordinates 
and momenta were generated by randomly sampling the 
phase-space distribution of the rotating Morse-oscillators 
and the spatial orientation of each molecule was taken 
from a uniform distribution within 4tt steradian. The 
magnitude of the angular velocity for a given vibrational 
phase was set to the angular momentum calculated from 
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Figure 7: Opacity functions for complex formation (green diamonds) 
and charge transfer (blue triangles), and probability curve for charge 
transfer (red triangles) if the complex is formed. Error bars corre- 
spond to a 2cr standard deviation. 



L = ( J( J+l)) 1 / 2 ^ and its direction was randomly sampled 
from a uniform distribution in the plane perpendicular to 
the axis of the molecule. 

Impact parameters b were between and 15 ao in steps 
of 1 ao, and in steps of 0.2 ao between 12 and 13 ao where 
the opacity function drops to zero (see below). The ini- 
tial separation of the center of mass of the two molecules 
was 20 A 38 ao), for which the intermolecular inter- 
action energy was less than 1.4 x 10 -4 eV. According to 
the Langevin-Gioumousis-Stevenson (LGS) model of ion- 
molecule reactions the rate does not depend on the rela- 
tive velocity of the partners [46]. Therefore, the relative 
velocity of the colliding molecules was set to the calculated 
beam flow rate (787 m/s) without dispersion. 

To validate the simulations, the probability of complex 
formation and the probability of CT were determined for 
a range of impact parameters b based on 500-1000 tra- 
jectories at each value (see Figure [7]). According to our 
quasi-classical model, the probability for CT was found to 
be 50% (within the statistical uncertainty) up to impact 
parameters b = 9 ao- The CT probability starts dropping 
slowly between 10 and 12 ao, and above 12 ao it decays 
steeply to zero around 13 ao- The probability of complex 
formation is 100% up to 9 ao and drops proportionally 
to the CT probability. Once the complex is formed, the 
probability for CT is « 50% regardless of the impact pa- 
rameter. This finding is expected based on symmetry ar- 
guments and considering that rapid charge-redistribution 
takes place before the complex decays. It is also in agree- 
ment with previous experimental results [25] . 

From the integrated opacity function an integral cross 
section for CT of a tot = (243 ± 19) was obtained. Mul- 
tiplication with the relative velocity (787 m s _1 ) yields 
k = (5.36 ± 0.42) x 10" 10 cm 3 s _1 for the second order 
rate coefficient for CT. This value compares favorably with 
previous experimental results of 4.24 x 10 -10 cm 3 s _1 , 
6.6 x 10" 10 cm 3 s -1 and 5.0 x 10" 
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Figure 8: Percentiles representing the distribution of lifetimes r for 
the Nj~ complex as a function of the impact parameter b (left panel) 
and integrated probability density p (right panel). The median life- 
time is around 10 ps, the complex decays within 100 ps in « 94% of 
all observed cases. 



Figure 9: Percentiles representing the rotational-state distributions 
for the N ^ products after breakup of the reaction complex as a func- 
tion of the impact parameter b (left panel) and integrated probabil- 
ity density p(iV + ) (right panel). The products are predominantly 
formed in rotationally excited states. 



for the CT reaction between 15 N+0 = 0) and 14 N 2 (> = 0) 

izasaasi. 

The lifetime of the complex was determined as the time 
difference between the last and the first crossing of the 
surface separation radius 7.09 ao (see Sect. 3.1 ). As Figure 
[^suggests, the median lifetime of N|~ is ~ 10 ps and does 
not or only slightly depends on b. Lifetimes up to 100 ps 
are found for 94% of the trajectories which qualitatively 
agree with previous estimates [24J [49] . 

4-4- Translation-to-rotation energy transfer 

The salient quantity which can be extracted from the 
trajectories to aid in the interpretation of the experimental 
results is the distribution of rotational angular momentum 
p(iV + ) of the N2" products after complex formation and 
decay. Figure [9] reports p(N^~) and suggests that for the 
majority (97%) of cases 7V + < 13. On energetic grounds, 
a total of 0.045 eV is available in collisional energy for 
the react ants, and the initial rotational quantum number 
of N2 can be as high as 3 (corresponding to a rotational 
energy of w 0.003 eV), which implies iV+ < 14. The N+ 
products are predominantly formed in excited rotational 
states, which supports the experimental findings of deple- 
tion of the ground state population caused by reactive 
collisions with the molecules in the beam. Note that the 
theoretical product state distribution cannot be compared 
directly with the experimental findings reported in Fig. [4j 
because the ions probed in the experiment result from 
a sequence of CT reactions, whereas the simulations only 
reflect single collision events. 

However, the trajectory calculations clearly reproduce 
the rotational excitation of the products observed in the 
experiment and reveal the mechanisms of the underlying 
T-R energy transfer occurring during CT. The reactant 
molecules approach each other in a random orientation 



and are accelerated towards the potential well of the lin- 
ear reaction complex. The torque exerted while forcing the 
reactants towards a linear configuration excites a counter- 
rotation of the N2 moieties which results in highly excited 
bending and torsional vibrations of the complex. More- 
over, anharmonic couplings lead to a practically complete 
redistribution of the available energy over all vibrational 
degrees of freedom during the lifetime of the complex. 
Upon its breakup, large- amplitude bending and torsional 
vibrations are converted into product rotations. Repre- 
sentative trajectories, projected on the PES, are shown in 



Fig. 10 to illustrate this effect. 



Close inspection of the trajectories reveals the likely 
presence of additional processes resulting in the rotational 
excitation of the products, including rotation-vibration 
coupling, back scattering and recrossing which all may con- 
tribute to the final state distribution. Under the present 
experimental conditions, the orbital angular momentum 
of complex-forming collisions is computed to be L < 140ft 
which is available for conversion into rotational motion of 
the reaction complex. Coriolis forces may lead to a cou- 
pling of the complex rotation to its internal motion, pro- 
viding another mechanism for the transfer of angular mo- 
mentum to the fragments. Furthermore, the trajectories 



shown in Fig. [10] suggest that a wide range of scenarios is 
possible. They include direct mechanisms, illustrated by 
the blue trajectory, and can range to bound states with 
almost full redistribution (randomization, i.e. IVR) of the 
internal energy as shown for the orange trace. This is rem- 
iniscent of the situation recently encountered in the vibra- 
tionally induced photodissociation of sulfuric acid where 
all regimes from prompt reaction to complete IVR were 
found, depending on the amount of internal energy made 
available to the molecule [50] . 
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R(a ) 



Figure 10: Projection of the PES onto the R and 62 coordinate (see 
Figure [6} with r\ = T2 = 2.08 ao,#i = 0, and ip = 0. Contour 
lines are drawn between and 1.262 eV. Three typical trajectories 
leading to rotational excitation of the products are reported. The 
arrows indicate the direction of motion (incoming and outgoing). 
The complex lifetimes amount to 87 fs (blue), 341 fs (green) and 
6604 fs (orange). 



5. Summary and conclusions 

In the present study we have demonstrated for the 
first time ion-molecule reaction studies with state-selected 
Coulomb-crystallized molecular ions using the symmetric 
CT reaction as an example. By rotational state 

selection of the reactant ions and their localization in space 
by sympathetic cooling in an ion trap, their internal and 
translational motions were completely controlled. By si- 
multaneously cooling the neutral co-reactants to the low- 
est rotational states and collimating their kinetic-energy 
distribution in a supersonic molecular beam, ion-molecule 
reactions were performed with an unprecedented degree of 
control over both the collision energy and internal states 
of the reaction partners. The analysis of the product-ion 
rotational state distribution yielded for the first time in- 
formation on T-R energy transfer occurring during CT. 
The experimental results were analyzed and interpreted by 
making contact with quasi-classical trajectory simulations 
on a full-dimensional PES. The simulations reproduced the 
experimental findings and yielded important insights at an 
atomistic level into the mechanisms underlying T-R energy 
transfer. 
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